Incommensuration Effects and Dynamics in Vortex Chains 
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We examine the motion of one-dimensional (ID) vortex matter embedded in a 2D vortex system with 
weak pinning using numerical simulations. We confirm the conjecture of Matsuda et al. (Science 
294, 2136 (2001)) that the onset of the temperature induced motion of the chain is due to an 
incommensuration effect of the chain with the periodic potential created by the bulk vortices. In 
addition, under an applied driving force we find a two stage depinning transition, where the initial 
depinning of the vortex chain occurs through soliton like pulses. When an ac drive is added to the 
dc drive, we observe phase locking of the moving vortex chain. 
PACS: 74.60.Ge, 74.60.Jg 



The static and dynamical behavior of elastic media in- 
teracting with periodic or random substrates is a topic 
of intense study as it relates to a wide variety of physical 
systems, such as friction where atoms move over a peri- 
odic potential created by the underlying stationary atoms 
0, sliding charge-density waves and vortex motion 
in superconductors with random quenched disorder or 
periodic hole or dot arrays . A simplified case occurs 
when the motion is confined in ID or quasi ID systems. 
An example of this is vortex motion in superconductors 
with nano fabricated arrays of channels |J|. There are 
also naturally occurring defects that can act as quasi-lD 
channels for vortex motion, such as grain ^0 1 1 1 and twin 
boundaries |]l2|. The vortices in fabricated channels ex- 
perience two types of pinning: the first from the intrinsic 
random disorder in the sample, and the second from the 
periodic potential created by the immobile vortex lattice 
surrounding the channels. Experiments and simulations 
in these systems have observed commensuration effects 
and interesting dynamical states ||,^ . 

Another intriguing example of quasi ID vortex matter 
is the vortex chain state which occurs in layered super- 
conductors with tilted magnetic fields as first imaged by 
BoUe et al. . Here a portion of the vortices align in a 
string with a smaller lattice spacing than that of the bulk 
vortices. Recently scanning Hall probe and Lorentz 
microscopy |l5| images have shown that the vortex chain 
states exhibit a remarkably rich variety of behaviors and 
structures for various regimes. The recent Lorentz mi- 
croscopy experiments found that for low temperatures, 
both the vortices along the chain state and the bulk vor- 
tices are stationary. For increasing temperatures the vor- 
tices along the chain disappear, which is speculated to 
be due to the onset of motion or oscillation of the vor- 
tex chains, so that the chain vortex images are blurred 
while the bulk vortices are still stationary It was 

conjectured in that this effect is due to an incom- 
mensuration effect where vortices in the chain move in a 
periodic potential created by the bulk vortices. Since the 
vortices in the chain have a different spacing than the pe- 
riodic potential, the chain is incommensurate and hence 



more weakly pinned than if the chain were commensu- 
rate. These images also showed that the blurring first 
occurs where there a mismatch accompanied by a topo- 
logical defect in the chain. Additional evidence for this 
scenario is that as the field is increased, the periodicity of 
the potential decreases and the incommensuration effects 
are enhanced, causing the disappearance of the chains to 
shift to lower temperatures, as observed jl^. Dodgson 
|l6| has proposed an alternative explanation in which 
the image disappears due the sudden formation of tilted 
chains. In addition to the intriguing physics of the vor- 
tex matter in confined geometries it is also of paramount 
importance to understand vortices in the chain state or 
along grain and twin boundaries, as these can act as weak 
lines where fiux motion occurs, leading to the breakdown 
of superconductivity in the chains prior to the bulk. 

In order to examine the conjecture that the disappear- 
ance of the vortex chain is due to an incommensuration 
effect, as well as to explore dynamical consequences of 
these incommensuration effects, we have performed nu- 
merical simulations of vortex chains coexisting with bulk 
vortices in the presence of weak quenched disorder and 
thermal noise. We also investigate the effects of applied 
dc and combined dc and ac external drives on the chain 
motion. Our results should be relevant for the vortex 
chain state where the vortices remain aligned in the z- 
direction so that much of the physics can be considered 
2D. In addition our results are relevant for vortex motion 
in artificial channels and grain or twin boundaries where 
there are additional vortices outside the channels. 

We consider a 2D sample with periodic boundary con- 
ditions in the x and y directions. We model the vortices 
as repulsive point particles which interact with weak ran- 
dom quenched disorder, applied driving forces, and ther- 
mal noise. A portion of the vortices are confined to move 
along a quasi-lD channel. The overdamped equation of 
motion for a vortex i is 



dr 

„ — f . — r^'" 



nrp 
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Here ry is the damping constant that is set equal to unity. 
The vortex-vortex interaction force is 
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FIG. 1. A portion of the sample with vortices (black dots) 
and trajectories (lines). The vortex chain is in the center, 
(a) T/Tm = 0.15. (b) T /Tm = 0.25, showing that motion is 
occurring along the chain but not in the bulk. 

^vv _ —J2'j'^i^iUv{rij), where the vortex interaction 
potential is Uy — ~Ayln{r). The force from random 
quenched disorder {[^ comes from randomly placed at- 
tractive parabolic traps of range rp with rp/a = 0.1, 
where a is the vortex lattice constant and the traps have a 
maximum force of /p = 1.0. The driving forces ff include 
a dc force f["-^ which is modeled as a uniform drive on 
the vortices in the a;-direction. We start from zero driving 
and with small increments increase the dc driving force. 
At each increment we average the vortex velocity Vx over 
a fixed time interval. An ac drive can also be applied in 
the form f/"-^ = Asin(a;t)x, where A is the ac amplitude 
and Lv is the frequency. The forces from the thermal noise 
come from random Langevin kicks with < f'^ (t) 
and < f'^{t)f'^{t') >= 2r]kBTS{t-t'). We measure tem- 
perature in units of Tm, the temperature at which the 
bulk vortices melt. For high- Tc materials at fields where 
the chain state is observed, Tm is around 70 — 80K and 
for the low- Tc materials with nano-channels, Tm would 
correspond to Tc. Lengths are measured in units of a, 
and forces in A^, the vortex- vortex interaction prefactor. 

We first consider the effects of temperature in systems 
with no external driving force to test the conjecture in 
We start the system at T/T,„ = 0.0 where the bulk 
vortices form a triangular lattice with lattice constant a. 
Additional vortices are added along the channel, giving a 
vortex spacing in the channel of a' < a. We monitor the 
vortex displacements for fixed time intervals of 10'' MD 
steps, and measure the vortex positions and trajectories. 
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FIG. 2. (a) Time series for the relative a;-position of a sin- 
gle vortex in the chain for T/Tm ~ 0.15 (light curve) and 
T/Tm = 0.25 (dark curve). For low T, the vortex is sta- 
tionary, while for higher T, the vortex jumps back and forth 
by about one chain lattice constant a. (b) The power spec- 
tra S{i') for the full time series of the vortex position. For 
T/Tm = 0.15 (lower light curve) the noise spectra is of low 
power and white. For T/Tm = 0.2 (upper dark curve) the 
noise power is much higher. The dashed line indicates a power 
law l/,r with Q = 1.6. 

For low T there is little motion of the vortices in the 
chain or the bulk as seen in Figure 1(a). In addition the 
incommensuration in the chain is stationary or pinned in 
a well defined location. For increasing temperatures there 
is a transition to a state in which the vortex chain exhibits 
a higher mobility than the bulk vortices, as seen in the 
vortex trajectory traces of Fig. 1(b) at T/Tm — 0.25. 

To elucidate the exact nature of the motion along the 
chain, in Fig. 2(a) we plot the position of a single vor- 
tex in the chain for the states in Fig. 1(a) and Fig. 1(b). 
For the stationary chain the vortex remains immobile 
with only small thermal fluctuations in position of mag- 
nitude less than a. For the mobile chain state the vortex 
does not move continuously, but instead moves in discrete 
jumps of magnitude nearly a. These jumps do not occur 
throughout the chain simultaneously, but run through 
the chain as a soliton like pulse with the vortices jump- 
ing sequentially. The discrete jumps occur due to the 
depinning of the incommensurations, which move ran- 
domly through the sample, causing each vortex to jump 
one lattice constant as the incommensuration passes over 
it. The discreteness of the jumps arises from the periodic 
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potential in which the chain vortices sit, created by the 
stationary ordered bulk vortices. Since the incommen- 
surations are much more weakly pinned than the bulk 
vortices, they depin at a much lower temperature. 

Figure 2(a) also shows that although there is motion of 
vortices along the chain, a single vortex does not diffuse 
but only moves back and forth by a single lattice con- 
stant. This can be understood by considering a single 
randomly moving incommensuration. When the incom- 
mensuration traveling in the positive a;-direction passes 
over a vortex, the vortex moves by a in the positive x- 
direction. The same incommensuration cannot, however, 
later cause the same vortex to move by an additional dis- 
tance in the positive a;-direction, but only in the negative 
x-direction when the incommensuration returns. Thus a 
single randomly moving incommensuration produces no 
vortex diffusion along the chain but rather an oscillation. 
If a series of incommensurations all moving in the same 
direction pass over a vortex, the vortex could move sev- 
eral lattice constants in one direction. This is prevented 
because there is no net shear in the system. The moving 
incommensuration sets up a strain field that prevents the 
incommensuration from traveling large distances in the 
same direction. As a result, the incommensuration does 
not diffuse freely through the chain but tends to drift 
back and forth about its original starting position. 

Using a Voronoi construction, we find that the incom- 
mensuration in the vortex chain coincides with a local 
defect composed of one 5-fold and two 7-fold coordinated 
vortices. It is the ID motion of this defect that occurs in 
Fig. 1(b). The images in Rcf. |l^] also indicate that the 
initial blurring in the chain occurs where there are local- 
ized topological defects in the vortex lattice. This agrees 
well with our finding that it is the motion of defects that 
is occurring along the chains. 

In Fig. 2(b) we show the power spectra for the mo- 
bile and immobile chain. For the immobile chain, the 
noise spectrum is white, while for the mobile chain, the 
low frequency noise power is several orders of magnitude 
higher, due to the correlated low frequency motions of 
the incommensuration, and there is now a 1//" signa- 
ture, with a = 1.6. For T > T„i, when all the vortices 
diffuse randomly, the noise spectra again becomes white. 

In simulations where the chain is commensurate with 
the bulk vortex lattice, such that a — a' , the chain and 
bulk vortex mobility are the same. Thus, the increased 
mobility of the vortices along the chain, as shown in 
Fig. 1(b), is due strictly to an incommensuration effect. 
When the vortex density along the chain is increased, 
lowering a' and increasing the number of incommensu- 
rations, the thermal depinning temperature of the chain 
decreases. The depinning and phase locking described 
below persist as the number of incommensurations in- 
creases. 

A further probe of the relative pinning strength of the 
incommensurations along the chain, compared to the 
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FIG. 3. (a) Velocity V vs applied dc drive fd for fd in 
the x-direction at T/Tm = 0.0. In the first portion of the 
curve, for fd < 1.0, vortices along the chain have depinned 
while the bulk vortices remain pinned, (b) The consecutive 
vortex position along the chain vs time for a fixed dc drive of 
fd/ fp ~ 0.15. Each curve corresponds to the trajectory of a 
single vortex in the chain. 

pinning of the bulk vortices, is the application of a dc 
drive. We consider the case with weak random pinning 
where, in the absence of a chain, there is a single well 
defined elastic depinning transition. In Fig. 3(a) the 
velocity-force curve for the incommensurate case shows 
two depinning transitions; the first is the depinning along 
the chain, and the second at 1.0 is the bulk vortex 
depinning. A similar two stage depinning has been ob- 
served in artificial channels and periodic pinning arrays 
Imaging experiments in the vortex chain state also 
find that the chain depins first under an external driv- 
ing force [|l5|. In Fig. 3(b) we plot the positions of the 
vortices in the chain as a function of time just above 
the depinning transition, showing that the motion is not 
continuous, but occurs in a pulse or stick-slip motion. 
Individual vortices remain stationary most of the time, 
and then jump one lattice constant in sequence. The dis- 
turbance or incommensuration moves much faster than 
the individual vortices. This is the same type of motion 
that occurs during the thermal depinning; however, in 
the driven case the incommensuration moves in only one 
direction. A similar type of ID soliton motion has been 
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FIG. 4. (a) Velocity V vs applied dc drive fd for 
T/Tm = 0.0 with an additional applied ac drive. The curves 
have been shifted up for presentation. Different ac amplitudes 
from bottom to top are A = 0.0. 0.02, 0.04, 0.06, 0.1, and 
0.12. (b) The width of the critical depinning force -qo vs A. 
(c) The width of the first step 771 vs A. 

seen for vortex motion in 2D periodic pinning arrays in 
simulations and experiments Q where the vortices 
effectively channel along the pinning rows. 

We have also investigated the effects of a combined 
dc and ac drive on the moving chain state. Since the 
motion along the chain occurs over an effective periodic 
substrate, interference effects should occur between the 
ac drive and the frequency induced by the dc drive. In 
this case, it is the incommensurations, rather than the 
individual vortices in the chain, that undergo phase lock- 
ing. In Fig. 4(a) we show a series of velocity-force curves 
for the moving chain for different ac amplitudes and fixed 
ac frequency. Here clear phase locking effects can be seen 
as steps in the velocity curves. The widths of the critical 
depinning force and the higher order steps vary nonmono- 
tonically with ac amplitude. If the phase locking steps 
are of the Shapiro step type, the nth step should vary as 
the Bessel function Jn{A) Q. In Fig. 4(b) and Fig. 4(c), 
the widths of the depinning force rjQ and the first step rji 
fit to Jo and Ji , respectively, showing that the phase lock- 
ing is best described as a Shapiro step similar to those 
seen in sliding CDW's Jg] as well as 2D vortex systems 
with periodic pinning [g| . The steps shown here occur be- 
low the bulk depinning transition, for fd/ fp < 1.0, where 
only the chain vortices are depinned. 

In conclusion, we have investigated vortices confined in 
ID embedded in an ordered 2D vortex lattice with weak 
pinning. The vortices in the chain are incommensurate 
with the periodic potential created by the bulk vortices. 
We confirm the conjecture in Ref. ||l^ that the smearing 
of the vortex chain state is due to an incommensuration 
effect, where for increasing temperature the incommen- 



surations thermally depin and become mobile before the 
bulk vortices do. The thermal motion along the chain oc- 
curs through a soliton like pulse which moves randomly 
through the sample. The vortices hop by one lattice con- 
stant in the positive or negative direction as the pulse 
moves through. The onset of motion along the chain 
can also be observed as an increase in the voltage noise 
with a 1//" characteristic spectra. Under application of 
an external dc drive, we observe a two stage depinning 
where the vortices in the chain depin first and form run- 
ning solitons. When an ac drive is added to the dc drive, 
we observe phase locking of the moving incommensura- 
tions in the form of Shapiro steps in the velocity-force 
curves. Our predictions can be tested through noise mea- 
surements in the vortex chain state under application of 
external drives or temperature, and also apply to vortex 
motion in artificial channels. 

Note: Just before submitting this work we became 
aware of recent experiments that find phase locking of 
dc and ac driven vortices in artificial channels |l^ . 

This work was supported by the US Department of 
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P. Kes and A. Tonomura for useful discussions. 
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